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Every lattice for which the bond percolation critical probability can be found exactly possesses a 
critical polynomial, with the root in [0, 1] providing the threshold. Recent work has demonstrated 
that critical polynomials can be defined on any periodic lattice. In general, the polynomial depends 
on the lattice and on its decomposition into identical finite subgraphs, but once these are speci- 
fied, the polynomial is essentially unique. On lattices for which the exact percolation threshold is 
unknown, the polynomials provide approximations for the critical probability with the estimates 
appearing to converge to the exact answer with increasing subgraph size. In this paper, I show how 
the critical polynomial can be viewed as a graph invariant like the chromatic and Tutte polyno- 
mials. Like these, the critical polynomial is computed on a finite graph and may be found using 
the deletion-contraction algorithm. This allows efficient calculation on a computer, and I present 
such results for the kagome lattice using subgraphs of up to 36 bonds. For one of these, I find the 
prediction pc = 0.52440572..., which differs from the numerical value, pc — 0.52440503(5), by only 
6.9 X 10"''. 



Percolation is the study of the formation of random 
clusters on lattices. Given an infinite lattice, we declare 
each bond to be open with probability p, and closed with 
probability 1 — p. The resulting random clusters grow 
in average size with p until we reach the critical point, 
Pc, at which an infinite cluster first appears. The de- 
termination of Pc is an unsolved problem in all but one 
dimension and on two-dimensional lattices that are self- 
dual 3-uniform hypergraphs pj . An example is shown in 
Figure [1^, where the shaded triangle can represent any 
configuration of bonds or sites. Critical thresholds on 
these lattices are given by the Ziff criterion jil, 



P{A,B,C)=PiAB,C) 



(1) 



where P{A, B,C) is the probability that all three ver- 
tices are connected through open paths, and P{A, B, C) 
is the probability that none are connected. Application 
of (HI to find a bond threshold results in a polynomial 
in p, with order at most equal to the number of bonds 
in the unit triangle. We may consider each bond of an 
n— bond triangle to have a different probability, resulting 
in a critical surface. 



/(pi, ...,Pn) = 



(2) 



where the function / is at most first order in any of its 
arguments. I refer to this property as "linearity" . Exam- 
ples of such critical surfaces include the square lattice. 



S{PI,P2) =1-Pl-P2, 

and the the honeycomb lattice, 



(3) 



Hipi,P2,P3) = ^ + P1P2P3 - P1P2 - P1P3 - P2P3 ■ (4) 

Note that critical surfaces may be multiplied by an ar- 
bitrary constant, or even another function with no root 
in [0,1], without changing the predicted bond thresh- 
old. Removing the former freedom by requiring the con- 
stant term to be +1, and the latter by demanding linear- 
ity, makes the critical polynomial unique. Consider now 




FIG. 1. A 3-uniform hypergraph. The shaded region need 
not be a simple triangle. 



the martini-A lattice, with the probability assignments 
shown in Figure [T^. If we set pi = 1 , this bond is con- 
tracted and its end vertices merged, the result being the 
honeycomb lattice with two bonds doubled in parallel. 
This doubled bond can be replaced by a single effective 
bond, so we get, 



i/(l-[l-p2][l-p3],P4,P5). 



(5) 



By setting pi — 0, we delete the bond, and the result is 
the square lattice with bonds doubled in series. Thus, we 
have. 



S{p2Pi,P3Pb)- 



(6) 



The A lattice is a 3-uniform hypergraph and its critical 
surface given by ([T]). The only way to reduce to the 
correct deleted and contracted surfaces and preserve the 
required linearity property is to set 



A{pi,P2,P3,Pi,Pb) = 

PlH{l - [1 ~P2][l ~P3],P4,P5) 
+ (1 -Pl)S{p2PA,P3Pb), 



(7) 
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which has the formal appearance of an average of the two 
special cases. This form is general for lattices that can be 
solved using ([T]) , and constitutes the deletion-contraction 
formula for critical surfaces of such systems. As another 
example, we find for the martini-B lattice (Figure [T|d), 

B{pi,P2,P3,P4) ^P4S{pi,l - [1 -P3][l -Pi]) 

+ (1 ~ P4)S{piP3,P4)- (8) 

The deletion-contraction algorithm can be used to extend 
the definition of the critical polynomial to lattices for 
which the critical threshold is not known exactly. For 
the kagome lattice of Figure [2^, we may write, 

K{pi,P2,P3,P4,P5,P6) = 

P4B{1 - [1 -P5][l -P6],Pl,P2,P3) 
+ {l-p4)A{pi,P2,P3,P5,P6), (9) 

where A and B are given by ^ and ([8]). Critical poly- 
nomials are found by setting all probabilities equal, and 
in this case K{p,p,p,p,p,p) — 0, i.e.. 



1 



3p^ - 6p^ 



12p'*-6/+p^ = (10) 

which gives the well-known 0, 0] approximation pc = 
0.52442971..., compared to the numerical 5| Pc = 
0.52440502(5). Thus, a critical polynomial may be de- 
fined on any lattice for which known lattices appear upon 
deletion and contraction of a bond, but there is no guar- 
antee that the resulting prediction will be exact. In fact, 
the definition extends to any lattice, for if either of the 
lattices appearing in the deletion-contraction formula has 
an unknown threshold, one simply applies the formula re- 
cursively and continues this procedure until known cases 
appear. This method is used to compute other graph 
invariants such as the chromatic and Tutte polynomials 

In the following, I call the subgraph of a lattice, on 
which different probabilities are assigned, the "base" of 
the process, which is then tiled to form the infinite lattice. 
For example, in Figures [5^, [^b. El IH are bases for the 
kagome lattice consisting of one, two, four, and six unit 
cells. Some known and conjectured properties of critical 
polynomials are as follows: 

1. They are unique (up to a trivial multiplication by 
a constant) once the base and tiling are specified 
d, In particular, they are independent of the 
bonds chosen to delete and contract at each step. 
In this sense, the critical polynomial is a property 
of a finite graph. 

2. If the single-cell prediction is not exact, it is gen- 
erally very close, usually within 10~^ of the nu- 
merically determined threshold [l3|- In addition, 
proper choices of larger bases lead to predictions 
closer to the exact answer, with accuracy increas- 
ing with base size. This is not proven, but the 
argument for it is the wide variety of examples in 
which it is clearly seen to be true. 



3. Conversely, if the critical polynomial gives the ex- 
act threshold for a base of a single unit cell, as is 
the case for self-dual 3— uniform hypergraphs, then 
any critical polynomial using a larger base makes 
the same prediction, i.e. the original polynomial 
always factors out. This is also a conjectured prop- 
erty but it would seem to be necessary for consis- 
tency, and I have found no counter-examples. 

4. In many cases, if the single-cell polynomial does 
not give the exact answer, then it can be shown Q 
that no critical polynomial derived from a finite- 
sized base will give the correct percolation thresh- 
old. This is true for the kagome lattice, as discussed 
below. 

Critical polynomials have been found for all the 
Archimedean lattices 

B E El- The polynomial ((TOl) 
was originally found by Wu using his "homogeneity" [3| 
assumption and recently, he extended the method to in- 
clude kagome subnets [Q!], with excellent results fTi|. In 
fact, the homogeneity approximation gives a prediction 
for the full g— state Potts critical frontier, but at present 
it appears limited to the kagome lattice and its subnets as 
it relies on a transformation from the triangular lattice 
to these kagome-type graphs. Although I focus on the 
kagome lattice in this paper, mostly because of the inter- 
est it has attracted over the years (e.g. I, [Hill), the 
method presented here may be applied on any periodic 
lattice for either site or bond percolation. 

In previous work, critical polynomials were found "by 
hand", with the (4,6,12) lattice and its 18— bond unit 
cell probably representing the limit of what one would be 
inclined to do this way (see the appendix of [1|). How- 
ever, the recursive nature of the algorithm makes it an 
ideal problem for a computer, and in this paper, I report 
the results of using a program to calculate critical poly- 
nomials on the kagome lattice for bases containing 12, 
24, and 36 bonds. At each step, the program chooses a 
bond, and finds the two graphs that result from its dele- 
tion and contraction. It knows a small number of graphs 
(e.g., triangular, honeycomb, and square lattices), and 
it repeats the deletion-contraction algorithm recursively 
until it recognizes all the lattices it has found. The output 
is a set of function definitions, like equations ([7]), ^ and 
dH) which can be evaluated in a computer algebra package 
to get the critical surface and then the critical polyno- 
mial. Of course, there are many issues to overcome in 
programming this scheme, and a more full account will 
be given elsewhere. Here, I show some examples of its 
operation on the kagome lattice. 



Base of 2 unit cells 

The first extension beyond a single unit cell base is 
that shown in Figure [21 in which we employ a base using 
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a) b) c) 

FIG. 2. a) a base for the kagome lattice using two unit cells, 
shown in two different shades. Each bond on the base has a 
different probability; b) a tiling for a); another tiling that is 
equivalent to b). 

two unit cells which are indicated by different colors in 
Figure [5^. The polynomial can be written in the factored 
form 

- (1 - 3/ - + 12/ - + /) 

X (-1 -p2 _ 2p^ + 10/ - 10/ + 3/). (11) 

We recognize the first term in brackets as the polynomial 
for the 6— bond base. The second term has no real roots, 
and thus the prediction here is the same as the one we 
found previously, pc — 0.52442971... . There is no other 
way to tile this base to give a different result. 

Bases of 4 unit cells 

Using the base consisting of 4 unit cells shown in Figure 
[3^, we may tile it in two different ways as indicated in 
Figures [SJj and Starting with Figure the critical 
polynomial can be written in the factored form 

- (1 - 3p2 - 6/ + 12/ - 6/ + /) X 
(-1 - p2 _ 2/ + 10/ - lOp'^ + 3p^) X 
(1 - 2p2 - 4/ + 7/ + 24p''^ - 28/ - 64/+ 
172/ - 184/ + 110p^° - 36p" + 5p^2) (12) 

Once again, the first term in brackets is just the 6— bond 
polynomial pO)) . and the others have no roots in [0, 1]. 
The prediction is again the same as the 6— bond estimate. 

Things finally become more interesting when we use 
the staggered embedding in Figure 131:;. In this case, the 
polynomial is 

1 - 6/ - 24p^ - 24p^ - 24p^ + 27p^ 
+ 552/ + 1056pi" - 1224p" - 8548^^2 

- 4872pi3 + 68568^^"* - 50664^^^ 

- 226650p^^ + 643944p" - 843684^^^ 
+ 684384p^^ - 368886p2° + 133152/^ 




a) b) c) 

FIG. 3. A 24— bond base for the kagome lattice with two 
inequivalent embeddings. 

- 31068/2 + 4248/3 - 259/^^ = 0, (13) 

and the solution on [0,1] is Pc — 0.52440672..., which 
differs from the numerical result by 1.7 x 10~^, a great 
improvement over the single-cell 6— bond case. There 
does not appear to be a way to get a better prediction 
using a 24— bond base. 

Bases of 6 unit cells 

There are many options for bases of 36 bonds, and I 
will discuss only three here. If we take the 6-unit cell ver- 
sion of the base shown in Figure [3^ with the embedding 
analogous to Figure [3)d, we once again find a polynomial 
in which (|10l) appears as a factor. Thus, for this base 
and embedding, we get the 6— bond estimate again. It is 
probably safe to assume that this trend continues with 
larger bases and embeddings of this type. 

Turning now to the base shown in Figure |4^, in which 
the embedding is indicated by the matching shapes on 
the external vertices, we get the polynomial, 

1 - 3/ - 12p5 - 20p^ - 60p^ - 132/ + 56p^ 
+ 684pi" -f 1440p" -f 2108p^2 -f 2052p" 

- 10452^1^^ - 68708]5i5 - 82980^1*^ + 280152p" 
+ 1316026pi*^ - 49980pi^ - 12878976/° 

+ 5124684/1 + 90816816/2 - 199458252/^ 

- 12979085/* -f 816398808/^ - 1939348056/^ 
+ 2677229528/^ - 2575935942/^ 

+ 1832168220/^ - 984362272/° 

+ 400507236/1 - 121897767/^ + 26954680/^ 

- 4096134/* + 382956/^ - 16617/^ = 0, (14) 

which has solution on [0, 1] Pc = 0.52440607..., differing 
from the numerical value by 1.1 x 10^^. A different base 
and embedding is shown in Figure|4}D and has polynomial, 

1 - 6/ - 2Ap^ - 14/ + 36p^ -f 39p^ - lOOp^ 

- 462pi° + 780p" 4583pi2 + 4812pi3 

- 9276p" - 71600pi^ - 85626^^*^ 

+ 312336pi^ + 1091146/^ - 509340^^^ 

- 9675936p2" + 5297340^^1 + 66607704^^2 

- 151097304/3 - 5319734/* + 610494828/^ 
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FIG. 4. 36— bond bases for the kagome lattices. Matching 
shapes on the external vertices indicate how each base is tiled 
to create the lattice. 



27 



1461237180^^" + 2022998000p 
1949295060p2^ + 1387593528/^ 
745850356p^° + 303533928p^^ 
92388675p^2 + 20427736p^^ - 3103578p^^ 
290052p35 



12579p^^ = 0, 



(15) 



with solution in [0, 1], Pc = 0.52440572..., slightly better 
than the previous estimate and within 6.9 x 10^^ of the 
numerical value. 

Clearly, we are observing the predictions converging to 
the exact answer with increasing base size. From these 
few examples it also appears that the estimates are con- 
verging from above, as they all seem to be greater than 
the numerical value. However, I know of no argument 
that guarantees this trend will continue. Note also that 
no polynomial for the kagome lattice derived from any 
finite-size base will give the exact answer. This can be 
seen as follows. Using the base in Figure H^, we delete 
many of the bonds to leave a single unit cell and a few 
connecting bonds, as in Figure [5^. Now, we may con- 
tract the connecting bonds to recover the kagome system 
in which the base is a single cell, as in Figure [5}d. By 
uniqueness, the prediction for this case must be equation 
([TU)) . which contradicts the prediction found by setting 
all bonds equal in the full critical surface, i.e. equation 
(fn)) . Thus, if the single-cell polynomial does not provide 
the exact threshold, and for the kagome lattice it does 
not, then, although we can get arbitrarily close by using 
ever larger bases, no finite critical polynomial derived in 
this way will ever solve the problem. 

I have presented critical polynomials for the kagome 
lattice, up to bases of 36 bonds. This is the limit of 
the current implementation of the algorithm, as it be- 
comes progressively more difficult to add bonds due to 
the exponential complexity. It is not uncommon for a 
large calculation to produce over a million function def- 
initions. Nevertheless, there is significant room for im- 
provement in the efficiency, as the rule used to choose 
the bond for deletion and contraction at each step can 



have a significant effect on the rate of reduction to the 
known simple cases. I have hardly explored this issue 
and presently use what amounts to a random bond se- 





b) 



FIG. 5. a) the subgraph and embedding resulting from delet- 
ing 22 bonds in Figure [4^, which is a single cell of kagome 
lattice and some connecting bonds; b) contracting the con- 
necting bonds gives the kagome lattice partitioned into unit 
cells. 



lection, rejecting only those choices that lead to undue 
complications. Moreover, this algorithm is perfect for a 
parallel implementation as it would require little inter- 
processor communication. It remains to be determined 
how much extra performance can be wrought from these 
considerations, but the problem will hopefully be seen as 
an interesting computational challenge. 

I thank Robert Ziff for helpful advice on this 
manuscript. This work was performed under the aus- 
pices of the U.S. Department of Energy by Lawrence Liv- 
ermore National Laboratory under Contract DE-AC52- 
07NA27344. 
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